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Using molecular dynamics simulation, we calculate fluctuations and response for steadily sheared 
hard spheres over a wide range of packing fractions <j> and shear strain rates 7, using two different 
methods to dissipate energy. To a good approximation, shear stress and density fluctuations are 
related to their associated response functions by a single effective temperature T c a that is equal to 
or larger than the kinetic temperature T^ n . We find a crossover in the relationship between the 
relaxation time r and the the nondimensionalized effective temperature T c g/pa 3 , where p is the 
pressure and a is the sphere diameter. In the solid response regime, the behavior at fixed packing 
fraction satisfies T7 oc exp(— cpa :i /T c ff), where c depends weakly on cj>, suggesting that the average 
local yield strain is controlled by the effective temperature in a way that is consistent with shear 
transformation zone theory. In the fluid response regime, the relaxation time depends on T e f[/pa 3 as 
it depends on Tkin/pcr 3 in equilibrium. This regime includes both near-equilibrium conditions where 
Teg ~ Tkin and far-from-equilibrium conditions where Toff 7^ Tki n . We discuss the implications of 
our results for systems with soft repulsive interactions. 



I. INTRODUCTION 

Many applications require understanding how disor- 
dered materials flow under an external load such as 
a shear stress [TJ [5]. Deriving such an understand- 
ing from thermodynamic principles requires identifying 
relationships between thermodynamic parameters and 
transport. Since external loads drive systems out of 
equilibrium, their thermodynamics cannot be described 
solely in terms of equilibrium parameters like tempera- 
ture and pressure. However, simulations [SHI] and exper- 
iments [1TJHT2] show that sheared spherical particles pos- 
sess an effective temperature that relates low-frequency 
fluctuations of various observable quantities to their as- 
sociated response functions. This effective temperature 
is likely to feature prominently in any thermodynamics- 
based description of flow in glassy materials [13]. First, 
the relationship between transport and effective temper- 
ature must be established. 

In Ref. pT] we begun to establish a relationship be- 
tween mobility and effective temperature T Q g by study- 
ing the behavior of a mixture of soft, repulsive disks in 
two dimensions at a packing fraction above random close 
packing over a range of shear strain rates 7 and kinetic 
temperatures T^ D . We found that the system only flows 
under shear if T g exceeds a threshold value similar to 
the dynamic glass transition temperature. Under strong 
shearing, we found that the average shear stress is pro- 
portional to a Boltzmann factor containing the effective 
temperature. These results imply that the effective tem- 
perature plays a key role in facilitating the mobility of 
sheared systems. 

Here, in order to further elucidate this role, we deter- 
mine the relationship between mobility and T e g over a 
wide range of parameters for a hard sphere model. The 
advantages of the hard sphere model are twofold. First, 
due to the hard-core nature of the interactions, T e g can 



only control the behavior of the model via a dimensionless 
parameter. Using the dimensionless parameter T c g/pa 3 , 
where p is the pressure and a is the sphere diameter, fa- 
cilitates a direct comparison with the equilibrium system, 
where T /pa 3 is the relevant control parameter. Second, 
the relationship between mobility and T e s/pa 3 should 
capture the leading order behavior of soft repulsive sys- 
tems, provided that parameters continue to be expressed 
in dimensionless form |15H17j . 

We find a crossover in the relationship between T c g and 
the relaxation time r that characterizes the mobility. In 
the solid response regime, r depends on T e g at fixed pack- 
ing fraction according to T7 oc exp(— cpa 3 /T c g ), where c 
depends weakly on <f>, suggesting that the average local 
yield strain is proportional to a Boltzmann factor con- 
taining the effective temperature. This is consistent with 
the effective temperature's role in shear transformation 
zone theory, where T c g controls the density of zones that 
are susceptible to shear deformation |13j . In the fluid re- 
sponse regime, the relaxation time depends on T c s/pcr 3 
as it depends on T k i n /p<r 3 in equilibrium. This regime in- 
cludes both near-equilibrium conditions where T e ff ~ T^m 
and far-from-equilibrium conditions where T e g 7^ T kin . 
This suggests that the mechanisms that control transport 
for the equilibrium fluid persist under shear even where 
the slow degrees of freedom responsible for transport fall 
out of equilibrium with fast degrees of freedom. 



II. MODEL AND METHODS 

Our model is a mixture of 4096 hard spheres of mass 
m, half each of diameter a and lAa to avoid crystal- 
lization. We conduct event-driven molecular dynamics 
simulations |18l I19j at fixed packing fraction <f> and uni- 
form shear strain rate 7 under Lees-Edwards boundary 
conditions. Maintaining such a system in steady state 
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requires dissipating energy. We perform two separate 
sets of simulations using two different methods to dissi- 
pate energy. For inelastic simulations, we dissipate en- 
ergy at each collision by imposing a coefficient of restitu- 
tion C < 1. For thermostatted simulations, we let colli- 
sions be purely elastic (C = f ) but dissipate energy into 
a thermal reservoir at temperature Tki n by periodically 
rescaling the velocities Hi to keep the kinetic temperature 
m(\vi\ 2 }i/3 within 1% of T^ n . Although there are three 
control parameters for each set of simulations ((f), 7, and 
C or TkinX one of them can be absorbed into the unit of 
time; for instance, time can be measured in units of 7 -1 
or cr-y/m/Tkin- There are therefore only two independent, 
dimensionless control parameters, the packing fraction <fi 
and a parameter that controls the strength of shearing. 
While 7<j-y/m/Tkin and C are both valid dimensionless 
control parameters, we choose to represent the strength 
of shearing in both sets of simulations by the ratio E/p 
of shear stress to pressure. For the thermostatted and 
inelastic simulations, E/p can be increased at fixed <f> 
by increasing the shear strain rate jcr^/m/T^ or de- 
creasing C, respectively; the equilibrium limit E/p — > 
corresponds to jaWrn/T^ — > and C — > 1. 

We characterize the mobility by the relaxation time 
t defined by Ar 2 (r) = cr/v3, where Ar z (i) = 

V (( r z W — r z(0)) 2 ) is the root-mean-squared displace- 
ment in the direction perpendicular to the shear plane. 

As we will discuss in the following section, we define the 
effective temperature T e g as the value that replaces the 
temperature T to satisfy two independent fluctuation- 
dissipation relations. The first fluctuation-dissipation re- 
lation is the compressibility equation relating the com- 
pressibility to the amplitude of density fluctuations, 



T XT = M = I Jim S (k). 

p p k^O 



(1) 



We define the isothermal compressibility out of equilib- 
rium as the compressibility at fixed Tki n and E/p, 



f 



XT 



dp 



(2) 



which we calculate by composing partial derivatives and 
taking finite differences. For a two-component fluid like 
ours the structure factor S(k) in Eq.[l]is the combination 
of partial structure factors [5D] 



S(k) = 



5 n (fc)5 22 (fc) - (5 12 (fc)) 2 
x\S 22 (k) + x 2 Sn(k) - 2y/xix 2 Si2(k) ' 



(3) 



where x a = 1/2 is the fraction of spheres of component 
a. We take the limit in Eq. [I] by calculating S(kz) for 
2irkL = 1,2,..., where z is the direction perpendicular 
to shearing and L is the periodic box length, and fitting 
the lowest-order expansion S(kz) — 5(0) + cfc 2 over the 
domain k < 0.3/<7. 



The second fluctuation-dissipation relation is the 
Einstein-Helfand relation 2 1 2.'i relating the shear vis- 
cosity r\ = E/7 to the amplitude of shear stress fluctua- 
tions, 



'/ 



1 ,. dH(t) 
— lim . 

T t^oo dt 



where 



H(t) 



2V 



(\G v (t) - G„(0) - EW| 2 > , 



(4) 



(5) 



G n = X^i m ^iJ/ii and V is the volume. The Einstein- 
Helfand relation is the analog of the Grccn-Kubo relation 



di(5E(t)££(0)) 



(6) 



for systems with discontinuous potentials. We calculate 
the limit in Eq. [4] by fitting the slope of H (t) vs t over the 
range 2r < t < 20r, where r is the relaxation time, so 
that we sample low- frequency fluctuations. We find that 
we must average over very long simulations to converge 
to the equilibrium expression, Eq. |4j for simulations with 
purely elastic collisions and no shearing. Using this con- 
vergence as a guide, we only compute the right side of 
Eq. [4] for simulations of duration at least 500r. 



III. SHEARED HARD SPHERES POSSESS AN 
EFFECTIVE TEMPERATURE 

Before analyzing the relationship between t and T c g, 
we first show that the effective temperature is a valid 
concept for sheared hard spheres by demonstrating that 
the two independent fluctuation-dissipation relations are 
satisfied by a single value of T e fj. 

We first consider the compressibility equation, Eq. [l] 
relating the isothermal compressibility to the amplitude 
of density fluctuations, quantified by the long-wavelength 
limit of the structure factor. In Fig.[l](a) and (b) we plot 
the left and right sides of Eq. [l] for both types of energy 
dissipation and five selected packing fractions as a func- 
tion of E/p. In plotting the left side of Eq. [l]in Fig. [l] 
(a), we replace T by the kinetic temperature Tki n - Two 
important features are apparent in Fig. [l] (a) and (b). 
First, except at very high E/p, the data for the inelas- 
tic and thermostatted simulations collapse. With one 
exception that we will discuss, we find such a collapse 
for all observable quantities that we measure, indicating 
that the behavior of the model is insensitive to the way 
that energy is dissipated. The breakdown of this col- 
lapse at very high E/p is due to an unphysical layering 
transition that occurs for sheared spheres coupled to a 
profile-biased thermostat [151 123 At these values 

of E/p, the data for the inelastic simulations should be 
considered the physically realistic branch. 

The second important feature of Fig. [I] (a) and (b) is 
that the data in panel (a) and (b) are not identical, in- 
dicating that the compressibility equation does not hold 
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FIG. 1: (a) Measured compressibility IkinXr vs shear stress 
E/p. (b) Compressibility expected from applying the com- 
pressibility equation (Eq. [Ij), S(0)/p, vs E/p. (c) Ratio of 
compressibility effective temperature to kinetic temperature, 
T x /Tki B , vs E/p, where T x is defined via Eq. [7] as the ratio of 
the compressibilities in panels (b) and (a). Data are presented 
for a selection of four packing fractions. Red, open symbols 
are for thermostatted simulations, while black, filled symbols 
are for inelastic simulations. 



under strong shear. In particular, while the measured 
compressibility at fixed (f> decreases with increasing E/p, 
the compressibility expected from applying the compress- 
ibility equation to the measured density fluctuations in- 
creases with increasing E/p. Nevertheless, following the 
procedure practiced in previous studies of effective tem- 
perature [3HH1 Ej , we can define a temperature-like pa- 
rameter T x as the quantity that replaces T to satisfy 



T _ S(0) 

X PXT ' 



(7) 



Fig. [T] (c) shows the ratio of this compressibility tem- 
perature to the kinetic temperature for the five pack- 
ing fractions as a function of E /p. For packing fractions 
4> = 0.55 and 0.58 below the colloidal glass transition [26] . 
we can easily conduct simulations in the near-equilibrium 
regime, where E/p is small and T x ~ T^ n . As E/p in- 
creases, T x becomes larger than T^ n , indicating that the 
degrees of freedom associated with compression are no 
longer in equilibrium with the degrees of freedom asso- 
ciated with velocity fluctuations. For packing fractions 
4> = 0.6, 0.62, and 0.63 above the colloidal glass transi- 
tion, all of our simulations are strongly nonequilibrium, 
with T x 3> Tkin- For each packing fraction, the left edge 
of the data represents the lowest strain rate jcr-^/ m/Xk in 
(or highest coefficient of restitution C) that we can ac- 
cess on the time scale of our simulations. Obtaining data 
for lower strain rates or higher restitutions of coefficient 
would require an intractably large separation between the 
time scales for collisions and shearing. 

In order for T x to be an effective temperature, it 
should satisfy more than one independent fluctuation- 
dissipation relation. The second fluctuation-dissipation 
relation that we consider is the Einstein-Helfand expres- 
sion, Eq. |4j relating the shear viscosity and the ampli- 
tude of shear stress fluctuations. We define the Einstein- 
Helfand viscosity as the right side of Eq. |4j replacing T 
with Tki n : 



1 



lim 

t— >oo 



dH(t) 
dt 



(8) 



In Fig. [2] (a) and (b) we plot 77 and ?7eh, respectively, as 
functions of E/p, again for both types of energy dissipa- 
tion and five selected packing fractions. We nondimen- 
sionalize the viscosities in Fig. [2] (a) and (b) by dividing 
them by a characteristic viscosity that is a product of 
the energy density p and the time scale y/m/pa; choos- 
ing such a characteristic energy density and time scale 
constructed from the pressure allows direct comparisons 
with systems with soft repulsions P"5"rfT?] . but our evalu- 
ation of the effective temperature is independent of this 
choice. 

Fig. [2] (a) shows a characteristic feature of the colloidal 
glass transition: for <f> = 0.55 and 0.58 below the colloidal 
glass transition, the viscosity is nearly uniform within a 
Newtonian regime at low E/p, while for cf> > 0.6 above 
the colloidal glass transition, the viscosity increases with 
decreasing E/p until it becomes too large to measure in 
our simulations. Outside of the Newtonian regime, the 
system exhibits shear thinning, with rjy/a/pm decreasing 
with increasing E/p. The apparent shear thickening for 
the thermostatted simulations at very large E/p is an 
artifact of the nonphysical layering transition. 

Even restricting the calculations to the state points 
for which we conducted simulations longer than 500r, 
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FIG. 2: (a) Measured shear viscosity r]^/a/pm vs shear stress 
E/p. (b) Einstein- Helfand shear viscosity t/eh \fa/prn vs E/p. 
(c) Ratio of viscosity effective temperature to kinetic temper- 
ature, T^/Tkin, vs E/p, where T,, is defined by the ratio of 
the viscosities in panels (a) and (b). Data are presented for 
a selection of four packing fractions. Red, open symbols are 
for thermostatted simulations, while black, filled symbols are 
for inelastic simulations. 



the data for ?7eh are much less precise. However, two 
observations are clear. First, the value of %h at given 
values of (ft and E/p is in most cases somewhat smaller 
for the thermostatted simulations than for the inelastic 
simulations. This is the only quantity that we measure 
that is sensitive to the energy dissipation mechanism at 
low or moderate E/p. We find that the discrepancy is ac- 
counted for by a somewhat faster relaxation of the shear 
stress correlation function in the thermostatted simula- 
tions compared to the inelastic simulations, though the 
time scales for the correlation function remain on the 




kin 



FIG. 3: Viscosity effective temperature vs compressibility ef- 
fective temperature for all simulations that were run long 
enough to measure both quantities. Red, open symbols are for 
thermostatted simulations, while black, filled symbols are for 
inelastic simulations. The solid lines representing T n — 3T X 
and T v = T x are visual guides. 



order of r for both inelastic and thermostatted simula- 
tions. Apparently, the profile-biased thermostat damps 
out some of the shear stress correlations that persist 
longer in the inelastic simulations. Second, analogous 
to our result for the compressibility, the behavior of ?7eh 
is markedly different than the behavior of rj at large E /p. 
While rj^a/pm uniformly decreases with increasing E/p, 
77eh \/ c /pm is non-monotonic for large <j>, first decreasing 
and then increasing with increasing E/p, and is roughly 
uniform for smaller (ft. 

We define the shear viscosity temperature T v as the 
parameter that replaces T to satisfy the Einstein- Helfand 
expression: 



kin 



r\ at 7?eh 



(9) 



Fig. [2] (c) shows that T n behaves qualitatively similarly 
to T x : for (ft = 0.55 and 0.58, T n ~ T kin at small E/p 
and increases at large E/p, while for (ft > 0.6, T v is sig- 
nificantly larger than Tk; n for all accessible values of E/p 
and increases uniformly with E/p. 

In Fig. |3j we show that T x and T n represent a consis- 
tent value of the effective temperature. We plot T^/Tkin 
vs T x /Tki n for a broad range of inelastic and thermostat- 
ted simulations within the range 0.2 < (ft < 0.636 and 
2 x 10~ 4 < E/p < 0.5. The points near T x /T kin = 
Tij/lkin = 1 represent near-equilibrium conditions, where 
the system is characterized by a single temperature, 
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kin 



T. For the thermostatted simu- 
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lations, the viscosity and compressibility temperatures 
away from equilibrium are strongly correlated, clustered 
around the line T n — T x representing a single effective 
temperature. We are unable to test this relationship 
for thermostatted simulations beyond T x /T^ n — 30, due 
to the layering transition at very large £/p. However, 
for the inelastic simulations, we are able to explore the 
range 1 < T x /Tki n < 300 without encountering artifacts. 
Throughout this range, we find a strong correlation be- 
tween T v and T x . However, the data clusters around a 
line representing T„ — 3T X , not T„ =T X . While we have 
already credited the discrepancy between thermostatted 
and inelastic simulations to the faster decay of shear 
stress correlations in the thermostatted simulations, we 
cannot explain the factor of 3 difference between T v and 
T x . However, given that the range (300) is much larger 
than the discrepancy (3), we find that T v and T x repre- 
sent a consistent value of the effective temperature to a 
reasonable level of approximation. 



IV. EFFECTIVE TEMPERATURE CONTROLS 
MOBILITY 

Having established that low-frequency density and 
shear stress fluctuations of sheared hard spheres are de- 
scribed by a consistent value of T G g , we now turn to estab- 
lishing the role of T e s in controlling the mobility, which 
we characterize by the relaxation time r. In the follow- 
ing, we focus on T x as the more precise measure of the 
effective temperature. Fig. [4] shows two different repre- 
sentations of the relationship between r and T x for the 
same five packing fractions and range of shear stresses 
examined in Figs. [T] and [2j In order to show the rela- 
tionship in a way that is independent of any arbitrary 
choice of units, we present the relationship in dimension- 
less form. Because there are no internal energy scales 
in the hard-sphere interaction, the effective temperature 
can only be compared the external energy scales like Tkin 
and per 3 . Since we are interested both in cases where 
T x ~ Tkin and where Tki n is negligible, we compare T x 
to per 3 rather than Tki n . In Fig. [i] (a) and (b), we repre- 
sent t using two different dimensionless groups. In Fig. g] 
(a), we plot the the average yield strain rj, the average 
amount of strain built up locally between relaxations. In 
Fig. [4] (b), we plot Ty/pa/m, the relaxation time nondi- 
mensionalized by the pressure. The time scale y/m/pa 
is the pressure-driven rearrangement time that sets the 
time scale for a loose configuration to be compressed into 
a close-packed configuration due to the pressure. Scaling 
both t and T e g by the pressure in Fig. [4] (b) serves to 
capture not only the behavior of the hard sphere system, 
but also the behavior of the related soft sphere system to 
leading order in pcr 3 /e, where e is the energy scale of the 
soft repulsion [151417) . 

In Fig. [4] (a) we examine the relationship between the 
effective temperature and the mechanical state of the sys- 
tem by plotting the average yield strain vs pa 3 /T x for our 
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FIG. 4: (a) Yield strain 7-7 vs pa 3 /T x . The lines are visual 
guides satisfying the equations rj — 10exp(— 0.22pa /T x ) 
and T7 = 1.8 exp(— O.llpa 3 /T x ). (b) Relaxation time nondi- 
mensionalized by the pressure Ty/pa/m vs pa 3 /T x . The solid 
curve is unsheared, equilibrium data for which T x = T. In 
each panel, data are presented for a selection of four packing 
fractions. Red, open symbols are for thermostatted simula- 
tions, while black, filled symbols are for inelastic simulations. 



selection of five packing fractions. We find a crossover 
between two regimes. In what we call the solid response 
regime, at low pa 3 /T x and high rj, the data at each 
packing fraction follows the relationship 



1-7 ~ c exp(-ciper 3 /T x ), 



(10) 



which appears as a linear relationship on the log-linear 
scale of Fig. [4] (a) . As illustrated by the visual guides in 
Fig. [4] (a), the coefficient Co ranges from 1.8 for <j) = 0.63 
to 10 for <f> = 0.55, and ci ranges from 0.11 for <j> = 0.63 
to 0.22 for cf> — 0.55. The appreciable average yield strain 



and its dependence on T x expressed by Eq. 10 indicates 
that within the solid-response regime, the system behaves 
like a continuously deformed solid with a mechanical re- 
sponse controlled by the internal state parameter T x /pa 3 . 
At high pa 3 /T x and low 7-7, the data for each packing 
fraction peel off of the linear relationship described by 
Eq. [lOj We will show shortly that these data correspond 
to a fluid response. The vertical parts of the data for 
= 0.55 and 0.58 indicate that for cf> below the col- 
loidal glass transition, the system approaches the equi- 
librium limit of 7 — > and pa 3 /T x — > pa 3 jT . We note 
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FIG. 5: Relaxation time r^/pa/rn vs (a) inverse compressibility temperature pa 3 /T x , (b) inverse kinetic temperature pa 3 /T^ n , 
and (c) shear stress E/p. Data are presented for a selection of four packing fractions. Red, open symbols are for thermostatted 
simulations, while black, filled symbols are for inelastic simulations. The solid curve that appears in panels (a) and (b) is 
unsheared, equilibrium data for which T x = T. In each panel, data are presented for a selection of four packing fractions. 



that the plot of E/p vs pa 3 /T x is qualitatively similar 
to Fig. [4] (b) , due to the fact that the shear modulus in 
a Maxwell-fluid model for hard spheres is roughly pro- 
portional to the pressure [T|>]. This is consistent with 
our result in Ref. [H] that E cx exp(— AE/T c g) for soft, 
repulsive disks above random close packing; the energy 
scale AE is controlled by the pressure. 

The form of Eq. [lO] is consistent with the shear trans- 
formation zone theory that describes plasticity of amor- 
phous solids [13] . In the theory, the effective tempera- 
ture describes the slow, configurational degrees of free- 
dom associated with rearrangements of neighboring par- 
ticles, and plastic rearrangements occur in zones whose 
internal configurational degrees of freedom are mod- 
eled as a system with two orientational states. The 
steady-state solution to the theory at low temperature 
is 7 = (l/2)eo-/Vzpz-R(E, T), where e is the strain 
per transformation, of order 1, Nz is the number of 
particles in a zone, pz — exp(— E'z/Teff) is the den- 
sity of zones, Ez is the energy to create a zone, and 
i?(E,T) is the rate of transformations in a zone [T5] , 
Without making any assumption about the form of 
i?(E,T), we assume that it is proportional to the dif- 
fusion coefficient; that is, we assume that the frequency 
of plastic rearrangements is linked to the frequency at 
which spheres make large displacements. This leads to 
7T = (l/2)Re Nz exp(— Ez/T e g), where R is the aver- 
age number of times a zone rearranges per relaxation 
time r. Comparing to Eq. 
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yields values of Re Nz be- 
tween 3.6 and 20 and values of Ez between 0.1 lpa 3 and 
0.21per 3 ; that is, the free volume associated with creating 
a shear transformation zone is approximately 0.1 la to 
0.22ct 3 . Our measured excitation energies Ez are simi- 
lar to the characteristic energy scales extracted from the 
dependence of Ty/pa/m on T/pa 3 for equilibrium hard 
spheres, for which fits of the Vogel-Fulcher and Elmatad- 
Chandler-Garrahan [371 HH] forms yield energy scales of 
0.18pa 3 and 0.25p<7 3 , respectively [T3] . 



In Fig. [4] (b) we demonstrate that for parameters out- 
side of the solid response regime, the relaxation time is 
controlled by the effective temperature in a characteristi- 
cally fluid way. In Fig. [2] (b) we plot t vs T x , where each 
quantity is nondimensionalized by the pressure, for the 
same set of data as in Fig. [4] (a). Fig. [4] (b) is a standard 
Arrhenius plot adapted to the effective temperature, a 
log- linear plot of TWpaJm versus the ratio of pressure 
to effective temperature, pa 3 /T x . For each packing frac- 
tion, the shear stress decreases to the right and upward; 
relaxation time increases and effective temperature de- 
creases as the shearing decreases. For comparison, the 
solid curve represents equilibrium simulations under no 
shear, which possess a single temperature T x = T^ n = T. 
In equilibrium, the relaxation time depends on the tem- 
perature according to the dimensionless function repre- 
sented by the solid curve, 



-VW^ = f(T/pa 3 ). 



(11) 



The packing fraction (f> increases as this equilibrium curve 
proceeds to the right and upward; the relaxation time in- 
creases as <j> approaches the colloidal glass transition and 
pa 3 /T increases according to the equation of state. We 
find that outside of the solid response regime-wherever 
the data in Fig. [4] (a) does not follow Eq. 10 -the data 
clusters around the equilibrium curve. Within this fluid 
response regime, the dependence of relaxation time on 
the two parameters <p an d E/p is captured by the equi- 
librium expression, Eq. 11 with T replaced by T x . This 



suggests that within the fluid response regime, the same 
mechanism that controls the relaxation of unsheared 
hard spheres also controls the relaxation of sheared hard 
spheres: spheres only flow if low-frequency fluctuations 
and large enough compared to the pressure to open up 
a sufficient amount of free volume |15j . However, for the 
sheared hard spheres, the low-frequency fluctuations are 
characterized by the effective temperature rather than 
the kinetic temperature. 
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EFFECTIVE TEMPERATURE AND THE 
JAMMING PHASE DIAGRAM 



FIG. 6: Color plots of (a) T x /T kin and (b) TyfpcT/ m as func- 
tions of kinetic temperature Tkin/po" 3 and shear stress S/p. 



In order to more clearly establish that the effective 
temperature, not the kinetic temperature, controls the 
mobility in the fluid response regime, we compare the 
dependence of Ty/pa/m on per 3 /T x with its dependence 
onpcr 3 /Tkin in Fig.[5](a) and (b), respectively. To aid the 
comparison, we plot pa 3 /T x and per 3 /T^ on the same 
logarithmic scale. Note that Fig. [5] (a) is a log-log version 
of the log-linear Fig.|4](b). For reference, Fig.[5](c) shows 
the same data plotted vs E/p. The equilibrium curve 
that appears in Fig. [5] (a) and (b) is identical, due to 
the fact that T x = T^ in in equilibrium. While much of 
the data-including large portions of the <j> — 0.55, 0.58, 
and 0.6 data and some of the (f> = 0.62 data-cluster near 
the equilibrium curve in Fig. [5] (a), the sheared data in 
Fig. [5] (b) only approach the equilibrium curve in the 
limit E/p — > for <\> = 0.55 and 0.58. For the packing 
fractions above the colloidal glass transition, where the 
near-equilibrium condition E/p << 1 is never met (see 
Fig. [5] (c)), the sheared data remain distinct from the 
equilibrium curve. 



We have shown that the ratio of effective temperature 
to pressure, T x /pa 3 , controls the mobility of sheared hard 
spheres according to two different mechanisms in two 
different regimes. In the solid response regime, the ef- 
fective temperature controls the mechanical state of the 
system, as evidenced by the dependence of the average 
yield strain. In the fluid response regime, the effective 
temperature facilitates relaxation by doing work against 
the pressure, analogous to how the temperature facili- 
tates relaxation in the unsheared system. It is instructive 
to organize these regimes within the jamming phase di- 
agram [TBI HH] that describes the relaxation of repulsive 
spheres as a function of temperature, packing fraction, 
and applied stress. Ref. [TB] recast the jamming phase 
diagram in dimensionless form in terms of the param- 
eters Tkin/pc 3 , E/p, and pc 3 /e, where e is the energy 
scale characterizing the repulsive interaction. For hard 
spheres, e = 00 and the jamming phase diagram reduces 
to the hard-sphere limit, pcr 3 /e —> 0. In this limit the 
relaxation depends on two parameters, an equilibrium 
parameter T^ n /pa 3 characterizing the strength of ther- 
mal fluctuations and a nonequilibrium parameter E/p 
characterizing the strength of shearing. 

Fig. [6] locates the solid response, fluid response, near 
equilibrium, and jammed regimes in the dimensionless 
jamming phase diagram. The three panels in Fig. [6] rep- 
resent the same sets of simulations, a large collection 
of all inelastic and thermostatted simulations for which 
E/p < 0.4 and T\^ n /pa 3 < 0.1. Each panel presents a 
color map of a different quantity. Fig. [6] (a) is a color 
map of the dimensionless relaxation time T\/pa/m. The 
region where data appears is the unjammed region of the 
phase diagram, where the equilibration and relaxation 
times are short enough that we can run simulations to 
calculate r. The red region indicates very large relax- 
ation times that are near the limits of our simulations. 
The dynamic jamming transition, defined as the locus 
of points for which Ty/pa/m equals some large but ar- 
bitrary number (like 10 4 ), separates the jammed region 
at small E/p and T^/pa 3 from the unjammed region at 
large E/p and/or large T^/pa 3 . 

Fig. [6] (b) shows that the unjammed region includes 
both the fluid response and the solid response regimes. 
Fig. [6] (b) is a color plot of the ratio of the effec- 
tive temperature T x to the the temperature T T = 
pa 3 f^ 1 (r yj pa / m) defined by inverting the equilibrium 
(zero shear) relationship, Eq, 
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between relaxation time 
and temperature. A value of T X /T T = 1 indicates that 
the relaxation time is predicted by Eq. [II] with T replaced 
by T x , corresponding to a collapse onto the equilibrium 
curve m Fig. [| (b) and Fig. [5] (a). Notice that a large 
portion of the jamming phase diagram is blue, demarcat- 
ing the fluid response regime where T X /T T is equal to or 
not much larger than 1 and the relaxation time is rea- 



sonably well predicted by Eq. 11 with T replaced by T, 



In particular, much of the dynamic jamming transition 
occurs within the fluid response regime. 

A comparison with Fig. [6] (c) emphasizes that fluid 
response does not imply near-equilibrium conditions. 
Fig. [6] is a color plot of the ratio of the effective tem- 
perature to the kinetic temperature. Values near 1, rep- 
resented by blue, indicate that the system is near equilib- 
rium, with fast and slow degrees of freedom characterized 
by a single temperature. Values much greater than I in- 
dicate that the system is strongly out of equilibrium, with 
an effective temperature much greater than the tempera- 
ture of the high frequency velocity fluctuations. While a 
large portion of the unjammed region is within the near- 
equilibrium regime, this regime is distinctly smaller than 
the fluid response regime of Fig. [6] (b) . 

Taken together, the three panels of Fig. [6] show that 
while the dynamic jamming transition at nonzero shear 
stress is largely controlled by a competition between low 
frequency fluctuations and pressure in the same way as 
at zero shear stress, these low-frequency fluctuations are 
characterized by an effective temperature that may be 
much larger than kinetic temperature. 

The jamming phase diagram provides a framework for 
applying our hard-sphere results to systems with soft in- 
teractions. Moving away from the hard-sphere limit, the 
relaxation of soft, repulsive spheres varies continuously 
with the parameter pa 3 /e that characterizes the softness 
of the potential [T]3 |TS] . Even far from the hard-sphere 
limit, where the relaxation time is not adequately ap- 
proximated by its form at pa 3 /e — 0, the dependence 
of T\Jpa jm on T/pa 3 (with no shear) can be mapped 
onto the hard-sphere behavior by computing an effec- 
tive hard-sphere diameter from the structure of the soft- 
sphere system [T7] . This suggests that while the numeric 
form of the relaxation changes far away from pa 3 /e = 0, 
the mechanisms remain the same. We therefore expect 
that for pa 3 ft > the effective temperature continues 
to facilitate relaxation via the same two mechanisms as 
in the hard sphere limit. While the boundaries between 
the two mechanisms may shift at elevated per 3 ft, we ex- 
pect that the solid response will continue to dominate at 
low Tkin/pa 3 and high T,/p while the fluid response dom- 
inates elsewhere. Analysis of our earlier simulations |14j 
of sheared two-dimensional soft disks in terms of dimen- 
sionless quantities is consistent with this picture. 

VI. CONCLUSIONS 

We have shown that the ratio of effective temperature 
to pressure, T c ff/pa 3 , controls the mobility of sheared 



hard spheres throughout the unjammed region of the 
jamming phase diagram spanned by packing fraction and 
applied shear stress. The effective temperature char- 
acterizes both the amplitude of static density fluctua- 
tions relative to the compressibility and the amplitude 
of low-frequency shear stress fluctuations relative to the 
shear viscosity. At high shear stress and low kinetic 
temperature, the effective temperature departs signifi- 
cantly from the kinetic temperature that characterizes 
high-frequency fluctuations. In this far-from-equilibrium 
regime, the effective temperature, not the kinetic tem- 
perature, controls the diffusive motion of the spheres. 

We find that the mechanism by which T c g /pa 3 controls 
the mobility depends on where the system lies in the jam- 
ming phase diagram, crossing over between two regimes. 
For very low values of the kinetic temperature relative 
to the pressure, the sheared hard spheres respond like a 
solid. They maintain an average local yield strain that 
depends on the value of T e g/pa 3 . The functional form 
of this dependence is consistent with the effective tem- 
perature's thermodynamic role in the solid-based shear 
transformation zone theory, in which the effective tem- 
perature controls the density of localized zones suscep- 
tible to plastic deformation. In contrast, for moderate 
or large values of the kinetic temperature, the sheared 
hard spheres respond like a fluid. The relaxation time 
depends on temperature in the same way as does the 
unsheared, equilibrium system, except that the relevant 
temperature is the effective temperature, not the kinetic 
temperature. This suggests that, just as the temperature 
facilitates flow in the unsheared system by doing work 
against the pressure, the effective temperature mobilizes 
the sheared spheres by working against the pressure to 
open up free volume. 

By expressing our results in relation to the dimension- 
less version of the jamming phase diagram, we have il- 
lustrated how these mechanisms likely apply for systems 
with soft repulsive potentials. It remains to be seen to 
what extent these mechanisms hold for systems with at- 
tractions. 
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